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d\ ' 

I We study the non-Gaussianity generated during multiple-field inflation. We provide an exact 

psj I expression for the bispectrum parameter Jnl which is valid beyond the slow-roll regime, valid for 

certain classes of inflationary models. We then study a new, exact multi-field inflationary model 
^ , considering a case where the bispectrum grows to observable values at the end of inflation. We 

I show that in this case the trispectrum is also large and may even provide the dominant signal 

^ ■ of non-Gaussianity. 

O ' 1 Introduction 

u: 

(— I [ The inflationary paradigm predicts an almost scale-invariant spectrum of nearly Gaussian perturbations, in 

Qh' excellent agreement with observations [1]. On the other hand, many different models of inflation fit present 

Q . day data very well, and further information is needed in order to test and exclude most of them. Non- 

^ [ Gaussianity, a signal of the non-linearity of the perturbations, has emerged in recent years as an important 

^ ■ discriminator between various inflationary models [2,3]. 

There are various known ways to generate large non-Gaussianity within the inflationary paradigm. 

CN ■ Models of inflation with non-canonical kinetic terms are characterized by an equilateral shape of non- 

! Gaussianity, generated on sub-Hubble scales [4]. Non-local inflation can also lead to large non-Gaussianity 

\^ ' with local shape [5]. These examples are non-generic, since inflaton perturbations are usually close to 

' Gaussian at Hubble exit [6, 7]. On the other hand, a local shape of non-Gaussianity can be generated 

, during the evolution of fluctuations on super-horizon scales H. Models have been proposed to produce a 

^ ■ large non-Gaussianity with this shape after inflation, such as inflationary models with an inhomogeneous 

' end [9, 10], modulated preheating or reheating [11], and the curvaton scenario [12]. However all of these 

. methods require the existence of additional, very light scalar fields, which play no role during inflation but 

j> I have a significant effect on generating the primordial curvature perturbation after inflation ends. Arguably, 

■ models which generate large, primordial non-Gaussianity during inflation are therefore more economical. 

^ ' Multi-field models of inflation, in which the curvature perturbation evolves at superhorizon scales, are 
able, at least in principle, to generate large non-Gaussian fluctuations during inflation. Most of the previous 
analyses of primordial non-Gaussianity in the multi-field case focused on separable Ansatze for the scalar 
field potential, either sum separable (for two-fields see [13], for an arbitrary number of fields see [14,15]) 
or product separable [16], and rely on the slow-roll approximation. Examples in which /tvl can grow large 
during slow-roll inflation for these classes of potentials were given [17], as well as the general requirements 
that the initial conditions of the background inflationary trajectory must satisfy. Numerical analyses have 
also been made, going beyond the assumption of separability and/or slow roll [18, 19], but no example of a 
model giving rise to a large non-Gaussianity at the end of inflation has been found in those cases. 
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^ Different possible shape dependences of the bispectrum (3-point function) are considered in [8]. 
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In this paper, we consider models of inflation driven by multiple scalar fields with canonical kinetic 
terms. We study the properties of non-Gaussianity, parameterised by the quantity fi\fL characterising the 
bispectrum, generated by the evolution of scalar perturbations on super Hubble scales. We develop a new 
method that provides, for certain classes of models, analytical expressions for /^rj^ valid in regimes beyond a 
slow-roll approximation. Our method is based on the first-order Hamilton-Jacobi formalism developed by 
Salopek and Bond [20], which allows us to express inflationary observables in the multi-field case, without 
having to focus on a slow-roll regime (see [21] for a similar application of this formalism to the single- field 
case) . 

The main idea on which our arguments are based is to focus on inflationary multi-field trajectories in 
which the Hubble rate, and not the potential, is separable. This approach allows us to analytically study 
cases where the non-Gaussian parameter /tvl becomes large, in regimes in which the slow-roll parameters 
are enhanced after Hubble exit, but before inflation ends. In particular, we provide analytical equations 
that express the non-linearity parameter /nl, in terms of quantities that generalise the slow-roll parameters 
during inflation. By means of the same method, we are also able to find new, exact multi-field homogeneous 
solutions, that satisfy our Ansatz, for which a large amount of non-Gaussianity can be generated before the 
end of inflation. This allows us to apply the general formulae we find on a concrete example, that shows how 
the fisiL parameter, and also both non-linearity parameters which characterise the trispectrum, can grow 
parametrically large as a function of quantities that characterise the model. In some cases the trispectrum 
through gi\fL can give the dominant signal of non-Gaussianity. 

The plan of this paper is as follows: first we briefly recap the 6N formalism, which allows us to compute 
the evolution of perturbations on large scales. In section [3] we present a method to calculate /nl exactly 
for some models, by using the Hamilton-Jacobi equations. In section [H we consider some explicit multi-field 
models of inflation for which we find exact solutions, and in one case we show that the non-Gaussianity can 
become observably large by the end of inflation. Finally we conclude in section \5\ 



2 Brief recap of 6N approach to evolution of perturbations 

In this paper we use the so-called 5N approach to study the evolution of curvature perturbations. This 
formalism has been found to be convenient for studying non-Gaussianity produced during inflation, and has 
been implemented in several previous studies. In this section we give a very brief review of the formulae 
necessary for the following discussion. 

The starting point is the expression for the number of e-foldings of expansion, given by 

^ = JJ H{t)dt, (1) 

evaluated from an initial spatially flat hypersurface to a flnal uniform density hypersurface. The perturbation 
in the number of e-foldings, 5N, is the difference between the curvature perturbations on the initial and flnal 
hyper surf aces. We take the initial time to be Hubble exit during inflation, denoted by t*. The flnal time tc 
corresponds to a uniform density hypersurface just before the end of inflation. The curvature perturbation 
is given by [22, 23] 

C = SN = Y^ Nj6yju + IY1 Nij6^iJ^.u + ■■■ , (2) 
/ IJ 

where Nj = dN/{dipl) and the index / runs over all of the flelds. The power spectrum and bispectrum are 
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defined (in Fourier space) by@ 



27r^ 



(CkiCk^Cka) = {27rf5\k^+k2+k3)B^{ki,k2,h) 



(3) 
(4) 



From this, we can define three quantities of key observational interest, respectively the spectral index, 
the tensor-to-scalar ratio and the non-linearity parameter 



/nl 



d log 
d log k 
Vt 

5 



/-a3 1^3 



B(^{ki,k2,k:i) 



Qkl + kl + kl 



(5) 
(6) 

(7) 



Here, "P* is the power spectrum of the scalar field fluctuations and Vt = 8P* = 8H^ / {Air'^ Mp) is the power 
spectrum of the tensor fluctuations. Both the spectra are calculated at the end of inflation and we will 
ignore any evolution after this time. As defined above, /nl is shape dependent, but it has been shown that 
the shape dependent part is much less than unity [6,7]. The ideal CMB experiment is only expected to reach 
a precision of /nl around unity [24], so we will calculate the shape independent part of /nl, denoted by /^^ 



in [13,16]. Whenever the non-Gaussianity is large, |/nl| ^ 1 as in the examples we are going to consider, 
we can associate /^^ ~ /nl- This {k independent) part of /nl and the spectral index can be calculated by 
the 5N formalism. 



V. 



n( — 1 



f(4) 
JNL 



-2e* + 



2 Eij^jNjiNi 
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(8) 

(9) 
(10) 



where = —H^./H'^, and H^, = H{t = ti,). In the next sections, we will find expressions for the previous 
quantities, valid also beyond a slow-roll approximation, for certain classes of multi-field infiationary models. 



3 The method 

In this section we present a method which allows us to obtain simple and compact expressions for the 
spectral index and the parameter /nl) i'^ terms of quantities associated with the homogeneous solution 
under consideration. The main advantage of this method is that we do not rely on any approximation in 
order to find our results. The final formulae then allow one to analytically study interesting regimes in field 
space, that cannot be probed within the slow-roll approximation, usually implemented when studying the 
evolution of perturbations during infiation. 

We start in section [3TT] by reformulating the equations of motion as first order Hamilton- Jacobi equations, 
following [20]. By implementing a suitable Ansatz, this formulation allows one to define two quantities that 
play an important role for studying the evolution of fiuctuations, that we start discussing from section [3^21 

''Note that different sign conventions of /jvl have been used. We use the same sign convention as WMAP [1], which is 
opposite to that used in [13, 16]. 
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3.1 The homogeneous equations 



We consider a system of two scalar fields (p and x with canonical kinetic terms, minimally coupled to gravity. 
The action, setting Mp^ = 8ttG = 1, where G is Newton's gravitational constant, is 



S 



-9 



We focus on homogeneous, classical trajectories, adopting a FRW Ansatz for the metric and flat three 
dimensional spatial slices: 



-dr + a{ty {dxj + dxl + dxf) . 



Defining the Hubble parameter as H = the equations of motion are 



(12) 



■ dH 



f dH 

V 50 



+ 



f dH 
dH 



X 



(13) 
(14) 



d(j) dx 
We concentrate on solutions satisfying the following Ansatz for the Hubble parameter H((p,x)- 

H{<P,x) = H^^\cp) + H^^\x) , (15) 

that is, we demand that it can be split as a sum □ of two pieces, each one depending on a sing le field % 
Notice that this is different with respect to the requirement of separable potentials usually done in the 
literature [13]. Indeed, given that the Hubble parameter appears squared in the left hand side of eq. p3|) . 
it is easy to see that the potential cannot be written as a sum of pieces, each depending on a single field. 

This separability Ansatz allows to define two quantities that characterise each homogeneous trajectory 
that satisfies eqs (I13p -()14p. The first is the number of e-foldings at time tc, that can be expressed as 



H{t) dt 



H 



(1) 



H 



(2) 



dx, 



(16) 



where we define H^^ = 
homogenous trajectory, 



^^^^^ and analogously for x- The second is a quantity C, conserved along each 

dx 



C 



H 



(1) 



+ 



H 



(2) ' 



(17) 



whose definition resembles an analogous quantity used in [13]. 

While the number of e-foldings parameterises the evolution along a given homogeneous trajectory, C 
characterises the motion orthogonal to it. We will see these quantities are useful in order to analyse the 
evolution of fluctuations at superhorizon scales, in a way similar to what is done in [13]. The main advantage 



^Hence we follow an analysis related to [13]. In an analogous way, one could also consider a situation in which the Hubble 
parameter splits as a product of single-field pieces [16] . 

® In principle, the Hubble parameter _ff is a function H{ip^''\t), t), with yj*^^ = and (^''^^ — x- However, by means of eq. 
HSj and (fTil) . we have 



dt 



4 OH 



3 at/jC") 



^(^(fe) dt 



p{^) —const 



2 dH 



dt^ 



= 



So that we have always H — H{ip^''\t)) without additional dependence on time [20] 
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of our Ansatz (llSp for the Hubble parameter is that the evolution of perturbations can be studied without 
having to rely on any approximation such as slow roll. The formulae that we are going to determine, then, 
can be applied also in regimes where slow-roll conditions are violated: this fact will be of crucial importance 
for the applications that we discuss in the second part of the paper. 



3.2 Evolution of perturbations 

As explained in section [21 we intend to follow the evolution of the curvature perturbation between an initially 
flat slice at time (corresponding to Hubble exit) up to a final uniform density hypersurface at time tc- To 
calculate the scalar spectral index and the /nl parameter, as in eqs. ^ and (llOp . we have to compute the 
first and second derivatives of the number of e-foldings along the initial values for the fields (p^, and x*- In 
order to do this, we use the properties of the functions A^(t*, tc) and C{4)i,, x*)i ^-^id we closely follow the 
arguments developed in [13] transposing them to our case of separable Hubble parameter. We can write 



dN 



1 

+ 2 



H 



(1) 



H 



(1) 



H 



(2) 



hP d(t>c hP dxc 

,0r. 'AC 



(2) 



hP dxc 



hP d^c 

iXc ,dlr. 



dx*- 



(18) 



We then need the derivatives of cpc Xc along the initial values of the fields. The following relations hold 

9C dC 



d(j)c 



dC 



, dXc 



dC 



+ 



+ 



dx* 
dC 

dx* 



dxi 
dx. 



(19) 
(20) 



since c^c and Xc are functions of the conserved quantity C. Also, we have 



^ = -J- ^ = J_ .21) 

By definition, the time t = tc corresponds to a hypersurface of constant energy density, so that H{(j)c, Xc) = 
const in a flat universe. This implies 

^H^l^ + ^H(^) = (22) 

dC dC '^^ ^ ' 

On the other hand, differentiating condition (jlTI) along C one gets 

#c 1 , dx. 1 ^23) 



H^ll dC ^(2) 



By plugging (f23l) inside ([22]) one can extract 
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This, together with ()17p . is enough to get 

d4>c d(j)c dC 
86^ dC dSi, 



H 



(1) 



+ 



H 



J 



and analogous formulae for the other derivatives of the fields. At this point it is convenient to define the 
following quantities: 



7^ 




H 



(2)' 



H 



H 



(2) 
XX 



H 



(24) 



(25) 



and (5 = 6'^+5^. Although their definition resembles the corresponding one for the usual slow roll parameters 
e and rj, they do not coincide with them when taking a slow roll limit 0. We assume that the quantities 6 
and 7 are much smaller than unity at Hubble exit t = t*, in order to have a reliable expansion for 5N as 
explained in section ([2]). During inflation, by definition, we have to ensure that the quantity en = —H/H^ 
is smaller than one. On the other hand, 



en 



1 



H 



(1) 



25. 



(26) 



implying that 5 must remain small during the infiationary period. The quantities 7 can however become 
large during inflation: this property will be very important in the following discussion. 

By means of these quantities, we can write 



dXc 



5, J)^ 



5} 



d(j)c 

dXc 
dx* 



51 
5i 



HcSl 
5c 



t f5^\^ 



5i 



(27) 



(28) 



We now have all the instruments to compute the derivatives of the number of e- foldings. The first derivatives 
result 



ON 



1 



Hi'^ + Zc 



dN 



1 K 



(2) 



Zr 



where 



dx* 2\/ 
HP5t-Hi''>5} 



(29) 



(30) 



^ It is nevertheless simple to work out the relation with the slow roll parameters, defined as 

2 



1 fW 



W 



V 



4,4, 



W 



00 



w 



and analogously for x- In a- slow-roll regime, one finds 
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This quantity controls the evolution of the curvature perturbation at superhorizon scales, as expected in 
multi- field inflationary models. 

The calculation of the second derivatives yields 

frf ^ + 1 

- + (31) 





1 


7X 




2 






1 






2 


2 5} 



d^N 1 dZc 1 dZc 



^* 2H^^J5: 

- Zr 1 dZr , , 

(32) 



(33) 



dx^ d<\>^ 2i7^ 2i7^ 

We now have all the necessary quantities for evaluating the expressions discussed in section [2j 
3.3 Curvature perturbation and non-Gaussianities during inflation 

We can use eqs. ([9]) and (jlOp to calculate the scalar spectral index during inflation, as well as the non-linear 
parameter /nl- In order to do this, it is convenient to introduce two new variables 

n = , v = , (34) 

with the property that u + u = 1. By means of these and of the results of the previous section, the spectral 
index for the curvature can be straightforwardly calculated from eq. ([9]), giving 

u{\- ^uj +v(l- 



nc-1 = -45,- 4^ "-^^ ^ (35) 

st ^ si 



Notice that, since it depends on the parameters u and v, it changes with time also beyond horizon exit, due 
to the correlation between entropy and adiabatic perturbations. On the other hand, if (5^ and 7* are small, 
the spectral index normally remains close to one during the entire inflationary period. 

We pass then to calculate /^li the parameter on which we are more interested in. We start by introducing 
the quantity 

A — 

where we define 

= (^S^j^ + 5'^7^) /5 . (37) 
One can then express the value for the fj^\^ by plugging the previously obtained results in formula (jlOj) , 



^-!&'^a-ti. (30) 



finding H 



5 ^,4, _ J('-|")+g('-g")+Kg-»)'-^ . (3„ 



Notice the formal similarity with equation (68) in [13]. The meaning of the various quantities, however, is different with 
respect to that paper. Various numerical factors are also different. 
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This is one of the main results of our paper. We reiterate that this formula is exact and is not based on 
a slow-roll expansion. Notice that the first two terms of eq. (j38|) are suppressed by 6t or and these 
quantities are less than unity since they are calculated at horizon exit. Only the last term proportional 
to A is not suppressed by quantities evaluated at horizon crossing. It depends on the quantity 7^^, that 
can become much bigger than one during infiation, rendering the size of A large enough to enhance fisfi to 
observable values. Our formalism allows to study this regime in a reliable, analytical fashion, contrary to 
an analysis relying on a slow-roll approximation. It is also possible to have a large /tvl while all slow-roll 
parameters are small, if the ratio of some of the parameters is very large and the background trajectory 
curves in an appropriate way [17] (see also [25,26]). We note in agreement with [13] that if one of the fields 
has reached a minimum so <^ = or x = at the end of infiation then ^ = and the non-Gaussianity at the 
end of inflation will be small. Therefore any model of infiation with a separable potential or separable Hubble 
factor with a large non-Gaussianity present at the end of inflation must have both fields still evolving, and 
therefore the presence of isocurvature modes. It would therefore also be interesting to study the evolution 
of the perturbations after infiation, to see if this leads to an important change in the observables [27]. 

In the next section, we apply our general results to various examples of infiationary trajectories, showing 
in concrete examples that our formulae can analytically probe regions in field space where the non-linear 
parameter /^^ is enhanced to large values. 

4 Exact solutions 

Our system of first order equations, discussed in section \37l\ also allows one to find new exact infiationary 
trajectories in the multi-field case, provided that the scalar potential satisfies some conditions. There are 
very few examples of multiple field exact solutions in the literature [28]. In this section, we present two 
classes of new exact solutions, with the aim to study the evolution of non-Gaussianities in each example. In 
the first example, section 14. H we analyse a quadratic potential, showing that in this case the non-gaussian 
parameter always remains small during the infiationary regime. In the second example, section 14.21 
we consider a combination of exponential potentials. In this case, non-Gaussianities can get enhanced to 
observable values, depending on the choice of the parameters in the potential. 

4.1 Quadratic potential: no enhancement of /nl 

We first present an example of exact configuration that does not lead to enhanced non-Gaussianities, since 
the quantity A defined in eq. (|36|) does not increase much during inflation. 
Consider a quadratic potential 



W = A + B (j) + C x + D (j)'^ + E + F (t)x 



(39) 



where the constants B, C and F are related to the other quantities by 




(40) 



(41) 



(42) 



The equations of motion can be solved choosing a Hubble parameter of the form 



H = Ho + Hi<p + H2X 



(43) 
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with 

Hi = l(^lA + E + Dy h! = ^, Hl = ^. (44) 

Setting for simplicity, and without loosing generality, = 0, the equations have the solution 

cl){t) = cP,-2Hit, x{t)=X.-2H2t, (45) 
a{t) = aoe(^°+^i'^*+^^^*)*-(^?+^')*'. (46) 

Notice that, although the solution is well defined for all t, at sufficiently large times the Hubble parameter 
becomes negative and the universe starts contracting. This point is, in any case, well beyond the end of the 
inflationary era, and we do not need to discuss it in our analysis. Inflation ends when the parameter en 
given by 

enit) = -H/H' = ^ + , = 2 5 (47) 

[Ho + H,<p. + H2X.-2{Hf + Hl)t]' 

becomes unity. Notice that the equations of motion, and eq. ()47p . allows one to rewrite the expression for 
the potential in a suggestive form 

W = 3H'^ -2 {Hi + Hi) = (l - ^) (48) 



ensuring that during inflation this quantity remains positive definite. 

In order to rely on the expansion for the number of e- foldings discussed in section [21 we impose that the 
quantities 5^''^ are small. Prom eq. ()47p . we learn that this is ensured as long as Hi, H2 <^ Hq. Notice that 
the 7 parameters in this example vanish identically: this simplifies the expression for the spectral index, 
that reduces to 



X 



nc = 1 - 45* - 4 — — - (49) 

so — 1 is proportional to terms suppressed by small (5v,'s parameters. Another consequence of having 
j<t>,x = 0, is the fact that the quantity A cannot increase much during inflation. So, following the discussion 
of the previous sections, we do not expect to find any enhancement of non-Gaussianities. Indeed, it is 
straightforward to get the following expressions for the parameters u and v: ti = 5^/(5*, v = 5j/(5*. A simple 
calculation then provides 



2 6t6i ^ {6}u-6tvy sf^x 



-fNL = 

= 2 6^, (50) 



that then remains small and constant during inflation. This is not surprising since, from eq. (j45|) . one learns 
that the inflationary trajectory follows a straight line in field space. To finish the discussion of this example, 
notice that spectral index and f^L are connected by the relation — 1 = —^/nl- 



4.2 Exponential potential: enhancement of /nl towards the end of inflation 

We discuss here an exact inflationary trajectory, for a potential based on a combination of exponentials 
that is instead able to produce large values of /nl at the end of inflation. Noticeably, a potential with a 
similar form may be motivated in string theory, in the context of Kahler moduli inflation [29] . We start by 
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presenting the inflationary trajectory, and showing how it can enhance the non-Unearity parameter /nl- We 
discuss the physical interpretation of our findings at the end of this section. 
The potential we consider is 

W{(l), x) = Uo (^1- Aie-'"^ + A2e-^'^'^ - Bie-^^ + B2e~^^^ + e~'''l>-Px^ ^ (51) 

where the parameters a and /3 satisfy 

2 3 6^2 a2 3 6^2 , . 

a = -TT , P = ^ • (52) 

The equations of motion discussed in section [3TT] admit the following solution for H: 

The overall factor to the potential Uq can be freely chosen so that the amplitude of the scalar power spectrum 
matches the observed amplitude of perturbations in the CMB. 

The solutions for the scalar fields read (from now on Hq = ^/Uq/S and we have again set = 0) 




1 



a 



In 



Ai Ho t 



The scale factor results: 

a{t) = 00(6" 



a Ai Hq t 



1 

2c? 



P^BiHot 



2^ ^Hot 



(54) 



(55) 



Notice that the solution becomes singular at late times, when the scale factor vanishes and the field values 
diverge. The singularity, on the other hand, occurs well after inflation ends, and we will not need to discuss 
it in our analysis. The quantities 6 and 7 read in this case 



6'f' = ^ 



tt2 -I a2 

4 



5^ = 



Hq^ (3 B^ -2/3x 

4 



Hp a^Ax 
IF 2 



Useful relations are 



T 
H 



-a 



Ho (3^Bi 
H 2 



Ho 



Notice that inflation occurs while 



en 



■yX 



+ 25^ = 2{^ 

a 



(56) 
(57) 

(58) 
(59) 



+ 21^ 



< 1. 



(60) 



In this example 7''^ and 7-^ can become much larger than unity at the end of inflation if a and/or [5 are much 
greater than one. This corresponds to a break down in slow roll, although en remains smaller than unity 
during the inflationary era. 
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We can assume that at the initial time t = t^,, the gamma parameters are smaher than a and /3, so that 
H-^ ~ Hq. Using the formahsm developed in section O we call 



2 2 ^ ~ 2 



2\ 6f Hr 



a Hq 



The number of e-folding at time tc reads 



2a 




+ In 




(5J H^\^ 1 / 1 
6! Hi) 



.(61) 



Focussing on the regime of large a and in which the first two terms in the previous expression are 
suppressed, the main contribution to the number of e- foldings comes from the terms that go like I/7*, as 
long as (5c or 6c are not too small (more on this later). 

We now proceed to calculate the quantities necessary for estimating the non-linear parameter /nl- One 
finds that 



1 6'P -6^ 

2 ^ 26 



¥ ^^/mx U6X 



a 



and 



7 



a 



The expression for A, evaluated at t = tc is given by 



A 



1 



6t6} 
61 



a 



6x + I3^6'i> 



OL\/6^^^^6t^ 



6t6} . 



(62) 



(63) 



(64) 



Although the quantities 6t and 6c are smaller than one during inflation, on the other hand the parameters 
a and (3 can be large; this allows the quantity A to become large during inflation, enhancing the value of 
•^NL' show this with an explicit example. 



A concrete example 

As a simple, concrete example, we choose a regime in which the parameters a and [5 are both large, let us 
say larger than some arbitrary quantity i? ^ 1 with no physical meaning by itself. We note from ()52p that 
we are therefore required to have < and B2 < 0. In this case the potential does not have a minimum, 
but inflation still ends through en growing larger than unity, and we can trust our results in this regime. 
The potential would need modifications which apply after inflation in order for reheating to take place and 
these modifications may provide a minimum for the potential. However a study of this issue is beyond the 
scope of this work. We parameterise the values of the quantities 6t and 6c at the end of infiation as 



J- 6- = ^ 
2m ' 2n 



such that 1 — 

m n 



1. 



(65) 
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Suppose that the quantities m and n are not too big. Namely, they satisfy the inequahty m, n <C i?^, such 
that the expression for the number of e- foldings, eq. (|61|). is dominated by the terms proportional to the 
inverse of 7^,. The same inequality also implies that //q — -f^* — -f^c- We choose |7f | = 4^^, I7* | = 45^, such 
that we can write 

Ntot = 60 ~ 10 (p + g) . (66) 



In order to have 60 e-foldings, we have to choose p + q = 6. This fixes the initial values for the fields 
X-k, to the values 



The conditions (|65p imply 



20paMi 



, m 
20p-,/^, 



20q 



n 



Since m and n are much smaller than E?, we expect p q. We hence have the simple relations 

1 .,.1 ^ 1 



5t 



6} 



11 - r. 



X 



From 



'toil-' ^^Hot 

dH]) and (p9]) for the tensor-to-scalar ratio, and from (j35]) for the spectral index, we find 

8 1 , 
< IQ-^ , 



tot 



1 



+ 7? 
- -0.04. 



tot 



and 

(67) 
(68) 

(69) 
(70) 

(71) 

(72) 
(73) 



So we have a completely negligible level of gravitational waves and a red spectral index in agreement with 
present day observations [1]. 

It is interesting to notice that for this choice of initial conditions we can write 



X 



~ — log 

a 

~ 75 log 



Ntot - N{t) + 
Ntot - N{t) + y 



m 




1 




+ 




2a\ 


a 


n 




1 




+ 




2f3\ 


P 



(74) 
(75) 



where we used eq. (I6ip . N{t) indicates the number of e-folds at time t, and Ntot corresponds to the total 
number of e-folds (that as we mentioned we take of order sixty). We know that the terms inside the square 
roots, in the argument of the logarithms appearing in the previous equations, are much smaller than one 
since m, n <^ E?. This implies that in a regime in which Ntot — N{t) ^ 1 we can write 



Px{t) = a 4>{t) -log 



(76) 



Since x is proportional to we have a straight trajectory in field space, implying that we cannot have 
conversion of entropy into adiabatic modes [30] , and consequently the curvature perturbation remains almost 
constant at superhorizon scales. In this regime, for the arguments we have seen in the previous sections, we 
cannot produce large non-Gaussianities in this context. However, towards the end of inflation, Ntot—N{t) < 1 



12 




and the trajectory acquires a sharp turn, since the terms inside the square roots in eqs. ([71]) and ([75]) start 
to give sizeable contributions. Precisely in this regime the curvature perturbation evolves at superhorizon 
scales, and parametrically large non-Gaussianities can be produced, as we are going to see now. 

We can calculate u and v in this case, finding (from now on we neglect corrections of order 1/R) 



m 



u 



1 In 

2 2 n + m 



' 2 2n + m 



2(5?. 



(77) 



It is simple to check that the curvature spectral index given in eq. Q remains close to one during the 
entire inflationary period. 
The parameter A becomes 

1 

A ~ 



V2{ 



3 



a 



+ 



n 



The leading contribution to the quantity f^) is 



u 



When evaluated at the end of inflation, using (j77p we find 



A. 



(78) 



6 



/nl ^ -2V2 



mn 



2 

an 



a 



+ 



n 



(79) 



We stress that this formula provide only the dominant contributions to Jnl- It is valid in the case laj, > 
> 1, and m,n < . 

As an explicit case, we take a = — 100,/3 = 20, m = 6,n = 6/5 and the initial conditions to satisfy 
(I74p and (I75p with Ntot = 60. Then, independently of the values of Ai and we find f^L — 58 from 
the simplified formula ([7^ . and /atl — 53 from the full formula ([55]) . See figures [Tj and [21 Notice that the 
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Figure 2: Left plot shows the trajectory considered for the parameters given after eq. (I79p superimposed 
on a contour plot of the potential. The square on the trajectory indicates a point along the trajectory one 
e-folding before inflation ends as (j) and x I'oll towards zero. This shows that the fields roll much more quickly 
during the final stage of inflation, and the trajectory curves at the end as discussed in the text after ([76|) . 
The right plot shows the potential for the same parameter values. Notice that inflation ends on the plateau 
long before the potential becomes negative. 



results for f^L differ from each other by quantities of order l/i? = l//3as expected by the approximations 
we have made. At the end of inflation, the parameters are 7*^ = 29, 7-*^ = —13 and this shows that the 
slow-roll approximation has broken down, see footnote [71 

Let us end with some considerations regarding our findings in this concrete example. Our expression 
for /nl in eq. (I79p depends on the initial values through p and q, and this dependence drops out at leading 
order in a. Hence f^i is independent of the number of e-foldings and scale independently. This is easy to 
understand, since we have seen that, in our example, non-Gaussianity can be produced only towards the 
very end of inflation, when the trajectory has a sharp turn in field space. 

Notice that our requirements of final values for the quantities dt and 5^ (both much bigger than l/R) 
imposes fine-tuning constraints on the initial conditions, since at leading order ml/R the values for 7^ and 
7^ must coincide (see eqs. (|671) - (f68]) and recall that p q). Importantly though this kind of fine-tuning 
does not force us to choose a very special field value which would be inconsistent with the 5N expansion, 
specifically we have checked that we satisfy S> and similarly for the x field [33]. For more discussion 
of this point see [26] . 

We stress that the enhancement of non-Gaussianity occurs towards the end of inflation. What happens 
just after inflation is a model dependent issue; we cannot address this question within the approximations 
used in this concrete model. For our choice of parameters we notice that Jml is still increasing at the end of 
inflation, see flguredl but for other choices it may start to decrease before inflation ends. We stress that in 
the regime where |a|, \f3\ ^> 1 the slow-roll parameters will necessarily become much greater than unity by 
the end of inflation, which may correspond to |7'^'*| ^ 1: this is one reason by which fi^L can become large. 
Our formalism is at least in principle suitable to study the evolution of non-Gaussianity after inflation, when 
the parameter en becomes larger than unity. We intend to address this issue in a future publication. 

^For an example of a potentially scale-dependent local /jv_l see [26], and for a scale dependent equilateral /jvl see for 
example [31]. Observational prospects were recently considered in [32]. 
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Trispectrum 



In this regime where a ^ [5 > R and m,n > \/E? it is also possible to give compact expressions for 
the trispectrum (4-point function) non-linearity parameters. Unlike for the bispectrum there are two non- 
linearity parameters which in principle are observationally distinguishable, since they multiply terms with 
different shape dependence. Using the delta- formalism they are given by [15,34,35] 

T.IJK NijNjkNjNk 25 ZijK N.ijkN^iNjNk 

TNL = ; -3 , gNL = 77 Z 73 • (^0) 



We already have all the ingredients to give the full expression for t^vl, eqs. (j29l) - (l33|l . and in principle it 
is a straightforward but lengthy calculation to differentiate these to get the third derivatives as required by 
gNL- The full results are extremely cumbersome, but it is not hard to estimate the size of these terms, in 
the regime we are currently studying. Indeed it is simple to see, from eqs. (I29p - (I33p . that at leading order 
in a the term dN/{d(j)i,) dominates at linear order in the delta- expansion, while at second order is the 
term d'^ N / (dcp'^) that dominates. With some more calculations, using the methods of Sec. 13.21 we find that 
the largest term at third order, at leading order in a and evaluated at the end of inflation, is given by 



f3m?rv' 



Now also using the fact that /? ^ 1 we can write the two terms for the trispectrum in terms of as 

/g \2 5^ 3n — m,/6 \^ 

TNL - [-fNL ) , -9NL - ^ J ■ (82) 

The result for ttvl is as expected for a model where a single field direction generates the primordial cur- 
vature perturbation [34] as it is a lower bound since in general tnl > (S/tvl/S)^ [36]. For a discussion of 
observational prospects for the trispectrum see [37]. We note that if m ~ 1 then from ()65p n ^ 1, so the 
trispectrum through g^L will give the dominant signal of non-Gaussianity through a large, negative gNL- 



5 Conclusions and Outlook 

We have considered models of inflation driven by multiple scalar fields with canonical kinetic terms. We 
studied the properties of non-Gaussianity, parameterised by the quantity Jnl characterising the bispectrum, 
generated by the evolution of scalar perturbations on super Hubble scales. We developed a new method, 
based on a first-order Hamilton-Jacobi formalism, that provides analytical expressions for /nl valid in 
regimes beyond a slow-roll approximation. The main idea on which our arguments were based is to focus 
on inflationary multi-field trajectories in which the Hubble rate, and not the potential, is separable. This 
approach allowed us to analytically study cases where the non-Gaussian parameter Jnl becomes large, in 
regimes in which the slow-roll parameters are enhanced well after Hubble exit, but before inflation ends. 

By means of the same method, we also found new, exact multi-field homogeneous inflationary trajectories, 
for which a large amount of non-Gaussianity can be generated before the end of inflation. This allowed us to 
apply our general formulae on a concrete example, showing how the Jnl parameter can grow parametrically 
large as a function of quantities that characterise a speciflc model. In the example we considered, the slow- 
roll approximation strongly breaks down shortly before the end of inflation: /nl grows to large values, and 
it also turns out to be scale independent. Furthermore the trispectrum is also large in this case, and for 
some parameter values and initial conditions the flrst signal of non-Gaussianity could come from a large 
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gNL- This is the first exphcit example of a break down in slow roll after Hubble exit giving rise to /at/, ^ 1 
at the end of inflation. 

Our method and results can be applied to other contexts. Other multi-field models of inflation, besides 
the ones we considered, can have inflationary trajectories that satisfy our Ansatz of separable Hubble 
potential. All these models can be studied by means of our methods, helping to analytically understand 
regimes in which /nl can be enhanced to observable values, and possibly finding new ways to enhance 
non-Gaussianities. In some cases, depending on the form of the potential and of the inflationary trajectory, 
the evolution of non-Gaussianity can be studied also after inflation ends, since our technique allows to go 
beyond a slow-roll approximation. For this reason it would be especially interesting to find an example of 
a potential with a minimum and large non-Gaussianity, so that we could address the question of how /tvl 
continues to evolve after the end of the inflationary period. Also, one can straightforwardly extend our 
method to study multi-field inflationary models with non-canonical kinetic terms, first studied in [38]. This 
can allow to further explore the consequences of having non-canonical kinetic terms for the evolution of 
scalar fluctuations, besides what has been done so far in the literature [4,39]. We hope to return to discuss 
these issues soon. 
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